FIXED-POINT FORMS OF THE PARALLEL SYMMETRIC 
SANDPILE MODEL 
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Abstract. This paper presents a generalization of the sandpile model, called 
the parallel symmetric sandpile model, which inherits the rule of the symmetric 
sandpile model and implements them in parallel. We prove that although the 
parallel model produces really less number of fixed points than that by the 
sequential model, the forms of fixed points of the two models are the same. 
Moreover, our proof is a constructive one, which gives a nearly shortest way 
to reach a given fixed point form. 



1. Introduction 

Sandpile model (SPM) was introduced by Bak, Tang and Weisenfeld [0 as a 
paradigm to describe the self-organized criticality (SOC) phenomenon in physics 
and has a variety of applications in physics, mathematics, economics, theoretical 
computer science § § |, |, @, ||, || g g]| • The simplest SPM model is 
that the system starts from a single column configuration (n), then at each step, 
one column gives one grain to its right neighbor if it has more than at least two 
grains comparing to its right neighbor. It was proved |n], [l2[ Q that this model 
converges to only one configuration at which the evolution rule can not be applied 
at any column (this configuration is called fixed point). Furthermore, all reachable 
configurations (which are obtained from the initial configuration (n) by applying 
several times of the evolution rule) are also well characterized and its configuration 
space is a lattice. The system SPM has been modified and generalized in several 
aspects to satisfy each particular purpose. In the context of chip firing games, 
cellular automata and informatics systems, the SPM model with parallel update 
scheme (i.e. at each step, all applicable rules are applied in parallel) received great 
attention |^|, ^, In ||], Durand-Lose showed that the transient time to reach 
a fixed point is linear in the total number of grains n when the parallel updated 
scheme is used, whereas it is 0(n 3 / 2 ) when the sequential one is used. 

To make it closer to the real physical phenomenon, Formenti et al. JTo| and 
Phan jl6|], generalized SPM so that grains are allowed to fall on both sides (left 
and right). This generalized model is called symmetric sandpile model and denoted 
by SSPM . The model has no unique fixed point any more. While Formenti et al. 
investigated the model by considering its configurations without caring its positions 
(that is, they identify all configurations which are up-to a translation on a line), 
Phan investigated the model in addition to its positions and showed the furthest 
position (comparing to the position at which the initial column is situated). The 
authors characterized reachable configurations in [flo| (resp. forms of reachable 
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configurations in |p~6|| ) starting from a single column configuration. Furthermore, 
they showed that the number of hxed-point forms of the model is exactly [y/n\ ■ 

In this paper, we study the SSPM model using a parallel update scheme. We 
denote this by PSSPM for further ease of reference. We stress that unlike SPM 
and PSPM, SSPM and PSSPM can have multiple fixed points when starting 
from a given initial configuration (as usual, we only consider single column con- 
figurations). It is also remarkable that in PSPM it is difficult to characterize all 
reachable configurations and so it is the same result in PSSPM. In this paper, we 
concentrate on fixed points of the parallel scheme PSSPM. It is straightforward 
that the set of fixed points of PSSPM is contained in the set of fixed points of 
SSPM . We show as main result in this paper that although this containment is 
proper in general, the containment of their forms of fixed points is an equality. 
Therefore, we can obtain all forms of fixed points of the sequential model by using 
the parallel update scheme with less than time. The proof is long and involved, 
and it has been divided into several subparts for better understanding. Indeed, the 
paper is structured as follows. Section 2 represents definitions and results about 
the characterizations of reachable configurations as well as the time of convergence 
for the three models: SPM, PSPM and SSPM. Section 3 contains the main 
result. First, we give a precise definition of PSSPM. Then, we introduce three 
procedures which are the building blocks for the proof of the main result (Theorem 
. Finally, in Section 4 we draw the conclusion and present some perspectives for 
future research on the subject. 

2. SANDPILE MODEL AND SOME GENERALIZATIONS 

In this section, we first give some basic definitions related to integer partitions 
and unimodal sequences. Then we represent the results of the sandpile model 
(SPM) and two generalizations (PSPM and SSPM) investigated in §, 
PHI- 

Definition 1. Let n and k be positive integers. Then 

(i) An integer partition is a non-increasing sequence of positive integers a = 
(oi, d2> • • • > a k)> moreover if a% + • • • + = n then a is called a partition 
ofn. 

(ii) A unimodal sequence of length k is a sequence ofk positive integers [a\, a^, ■ ■ ■ 
such that there exists an index 1 < i < k satisfying the condition a± < a2 < 

• • • < Oi > dj+i > • • • > ak-i > a*. The quantities defined by 

k 

h(a) — maxjai}*^]^ and w(a) = ^""^ dj 

i=l 

are respectively called the height and the weight of a. 

(iii) The reserve of a sequence a — (a\,a2, ■ ■ ■ >cife); denoted by a -1 , is the se- 
quence (a^, afc_i, . . . , a\). 

iv) A nth power of a sequence of positive integer a, denoted by a n , is the 
sequence obtained by concatenating n times a. 

Given a unimodal sequence a and an index 1 < i < k, we denote 

a<i = (ffli, • • • , Oi-i) and a>i = (a,+i, . . . ,Ofc), 

a<i = (a\, . . . , ai-i,ai) and a>, = (a*, a^+i, . . . , afc), 
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and call them the strict left sequenceand the strict right sequence of a by i, the left 
sequence and the right sequence of a by i, respectively. 

A discrete dynamical system is described by its configurations and its evolution 
rule. A configuration b is reachable from another configuration a if & is obtained 
from a by applying several times the evolution rule, and we write a — > b. We usually 
consider the system starting from one configuration, called the initial configuration, 
and then we investigate the set of all configurations reachable from this initial one, 
and we call this set the configuration space of the system. By this way, the system 
is well defined by its evolution rule and its initial configuration. A fixed point (or 
stable configuration) of the model is a configuration reachable from the initial one 
and on which the evolution rule can not be applied. 

The sandpile system is a discrete dynamical system used to describe the self- 
organized criticality (SOC) phenomena in physics. Our studying models are in 
a subclass of the sandpile system. In these models, a configuration, also called 
sandpile, is represented by an integer sequence (a\, 02, . . . , fflfc), the part ai is called 
the height of the pile (or column) i. In this paper, we always assume that for all 
sandpile models, the initial configuration is a single column containing n grains and 
the position at which these grains is situated is called its initial column. Therefore, 
each model is well-defined by its evolution rule. Wc first give here the definition of 
the simplest sandpile model (SPM) introduced in |J. 

Definition 2 (||). The Sand pile model is a system defined by the following SPM 
rule (right rule): 

• column i is right collapsible if ai — aj+i > 2 and when it collapses on the 
right it gives one grain to its right neighbor; 

• at each step, there is at most one collapse. 



Figure 1. SPM rule 



We denote by SPM(n) (resp. SPM) the configuration space of the sandpile 
model starting from (n) (resp. from any single column configuration). Reachable 
configurations of SPM are characterized as the following: 

Theorem 1 ([fj"lf). Let c be an integer partition. Then c is a configuration of SPM 
if only if it does not contain any subsequence of the form (p,p,p) or (p,p,p — l,p — 
2, q + l,q, q) for positive integers p, q satisfying < q < p. 

From this characterization, the fixed point of SPM is given by an explicit for- 
mula. 

Corollary 1 (|TT|). Given a positive integer n. Then SPM(n) has a unique fixed 
point which is of the form (p,p — l,p — 2, q,q,q — 1, ....2, 1) with p, q € N and 
q < p. Furthermore, the time to reach this fixed point is 0(n 3 ^ 2 ). 
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Figure 2. Configuration spaces SPM (6) and SPM (30) 

Now, we define the SPM in parallel (so called parallel sandpile model) introduced 
by Durand-Lose [§). 

Definition 3 (||). TTie parallel sandpile model is a system defined by the following 
PSPM rule: 

• at each step, each column collapses at most once; 

• at each step, all columns which are right collapsible collapse on the right. 

We also denote by PSPM(n) (resp. PSPM) the configuration space of the 
parallel sandpile model starting from the single column configuration (n) (resp. 
any single column configuration). 
Remark: 

- In SPM and PSPM, the first column is always a highest column. 

- The SPM is non-deterministic (since at each step, there may exist several 
collapsible columns) whereas the PSPM is deterministic. Both models 
have the same unique fixed point. 

- Although it is impossible to reach all reachable configurations of the se- 
quential SPM by parallelism, it is possible to get its fixed point in linear 
time by parallelism. 

The following result gives us the time to reach the fixed point by parallel model. It 
is 0(n) comparing to 0(n 3 / 2 ) in sequential model. 

Theorem 2 (§). Let n be a positive integer. The time to reach the fixed point 
from (n) in PSPM(n) is 0(n). 

Next, we represent a recent generalization of SPM where columns not only 
collapse on the right but also on the left. 

Definition 4 ( (To], |l6|). The symmetric sandpile model is a system defined by 
SSPM rule as follows: 

• addition to the right rule in SPM model, there is also the left rule, that 
mean one column i can give one grain to its left neighbor if ai — a^-i > 2. 

• at each step, there is at most one collapse. 
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(a) Configuration (b) Configuration space 

space 5PM(6) PSPM(6) 



Figure 3. 

We denote by SSPM(n) (resp. SSPM) the configuration space of the symmetric 
sandpile model starting from (n) (resp. any single column configuration). 
Remark: 

- The SSPM is a non-deterministic model since there may have columns 
which are collapsible on both sides and it may have more than one fixed 
point. 

- Unlike SPM, each configuration of SSPM is a unimodal sequence. 

- Unlike SPM and PSPM where the initial column is always a highest one, 
the position of the highest column of SSPM can be changed during the 
evolution. 




Figure 4. Configuration space SSPM(5) 



Figure [| illustrates all the reachable configurations of SSPM starting with 5 
grains, where the initial column is in the parentheses. Then, SSPM(5) has 4 fixed 
points as shown in the figure. It is also important to consider the fixed-point forms 
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which does not mention to their positions. Hence, we identify all fixed points which 
are up-to a translation. For instance, 11(2)1 and 1(1)21 have the same form 1121 
(but different in the position); 12(1)1 and 1(2)11 have the same form 1211. Hence, 
fixed points of SSPM (5) are in one of two forms 1121 and 1211. The study of form 
of SSPM is investigated in [10 and that of position is in Jl6| . In fact, these papers 



give a characterization of the form of SSPM. 

Theorem 3 (0, 0). A unimodal sequence a = (oi, 02, Ofc_i, Ofc) is the form 
of a reachable configuration of SSPM if only if there exists an index 1 < i < k 
such that the sequences a>i and (a<») are configurations of SPM. 

Furthermore, the enumeration of fixed-point forms of SSPM is also given. 

Theorem 4 ([[To|). The number of fixed-point forms of SSPM (n) is [\/n]. More- 
over, if P is a fixed point of SSPM(n) then P is of height either [\/v\ or [\/n\ — 1. 

3. The parallel symmetric sand pile model 

In this section we introduce another generalization of the sandpile model. In this 
model, we inherit the rule of the symmetric sandpile model and implement them 
in parallel. 

First we give precisely its definition. Like the other generalizations of the sandpile 
model represented in the previous section, we always start with the single-column 
configuration. 

Definition 5. The parallel symmetric sandpile model is a system defined by the 
following PSSPM rule: 

• At each step, all collapsible columns collapse; 

• For columns which are collapsible on both sides, it must choose exactly one 
direction to collapse. 




Figure 5. Configuration space PS SPM (5) 

We denote by PSSPM(n) (resp. PSSPM) the configuration space of the par- 
allel symmetric sandpile model starting with (n) (resp. any single column configu- 
ration). 
Remark: 
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(1) Unlike PSPM which is deterministic, the PSSPM is non-deterministic 
since although columns collapse at the same time at each step, there may 
have two directions (must choose one) for one column collapsing; 

(2) Since an evolution step by PSSPM rule can be considered as a combination 
of some evolution steps by SSPM rule, each configuration of PSSPM is 
a SSPM configuration. Furthermore, the configuration space of PSSPM 
is a subspace of that of SSPM and the set of fixed points of PSSPM(n) 
is a subset of that of SSPM(n). 

We notice that the set of fixed points of PSSPM is a proper subset of that 
of SSPM. Actually, SSPM(5) has 4 fixed points: 12(1)1, 11(2)1, 1(2)11, 1(1)21 
(see Figure |), but PSSPM(5) has only 2 fixed points: 1(2)11 and 11(2)1 (see 
Figure ||). However, one can observe that SSPM (5) has only 2 fixed-point forms 
as PSSPM(5), which raises a question about the correlation of fixed-point forms 
of the two models. 

The main result of this paper is to state that the set of fixed-point forms of 
PSSPM and that of SSPM are the same. Moreover, we can show an explicit 
evolution by PSSPM rule to reach any given fixed-point form of SSPM. 

Theorem 5. The set of fixed-point forms of PSSPM (n) is equal to that of SSPM (n) 
Consequently, there is [\/n] fixed-point forms of PSSPM (n) . 

As our proof lengthens in many steps, we break the proof into some lemmas. 
First, we sketch the main idea of the proof, we leave its details at the end of this 
paper after presenting necessary procedures. 

Main idea. For a fixed point P of SSPM(n), we construct a sequence of PSSPM 
transitions to obtain P from the initial configuration (n) . Because we are interested 
in the form of P but not in its position, we can suppose that the center column of 
P is at position (the notion of "center column" , one of highest columns, will be 
given later). In the constructed evolution, the column is always a highest one, so 
the choice of PSSPM rule in each step is in fact the choice of direction collapsing 
of the column 0. 

For a symmetric fixed point P, i.e. (P<o) 1 = P>o- The evolution is an Al- 
ternating Procedure, described as follows: at odd steps, the column collapses on 
the right, and at even steps, it collapses on the left. From (n) this procedure will 
converge to the symmetric fixed point P (see Corollary ||) . 

For P not symmetric, we can suppose that the column is the center of P, i.e. 

d = \w(P >0 ) - w(P <0 )\ = min| W (P >l ) - w(P <t )\. 

i 

Not loosing generality we assume that w(P > o) — w(P < o) > and let d — w(P > o) — 
w(P < o). The evolution by PSSPM rule is composed of three procedures: 

i) Pseudo-Alternating Procedure: a procedure from (n) to the configuration 
Q = (l,2,...,d- 1, (n - d 2 ),d,d - l,...,2,l) (see Lemma §). Note that 
w iQ>o) ~~ w (Q<o) is exactly d. 

ii) Alternating Procedure: a procedure from Q to the configuration R on which 
we could not apply any more the Alternating Procedure. It will be proved 
that R is of height h (see Lemma 0) . 

iii) Deterministic procedure: a deterministic procedure from R to P, where at 
each its step there is no any column collapsible on both sides (see Lemma 
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Like the symmetric case, we claim that at the end of the evolution we obtain P by 
Lemma ^, Lemma ^ and Lemma ||. □ 

To present the proof we hrst introduce some necessary procedures. The following 
procedure is implemented on integer partitions. 

Definition 6 (Atom Procedure). Let t be a non-negative integer and let a be an 
integer partition. The Atom Procedure of t steps is a sequence of t transitions 
starting from a described as follows: 

(i) The PSPM rule is applied at all steps. 

(ii) At each odd step one grain is added to the first column and at each even 
step no grain is added. 




FIGURE 6. 9 first steps of Atom Procedure from (4, 3,2, 1). The 
dark grains are added at odd steps 

We denote by Atom 1 (a) the configuration obtained from a after t steps of Atom 
Procedure. 

Recall that each sequential sandpile model is coded by a chip firing game on 
line by taking differences between two consecutive parts of each integer partition. 
By this we associate each integer partition a — (oi, . . . , a&) to the sequence of its 
differences d(a) = (di, . . . , <2&), where di = ctj — a^+x (by convention ak+i = 0). 
Therefore if a* — cn+i > 2 then di > 2. Moreover, if a' is obtained from a by one 
transition by SPM rule at position i then d(a') = (di, . . . , cij-i + 1, di — 2, di+\ + 
1, . . .,d k ). 

Denote by s(k) = (k, k — 1, . . . , 2, 1) the stair of k steps. The following lemma 
describes all configurations which are reachable from s(k) by Atom Procedure. 

Lemma 1. Let k,t be positive integers such that < t < 2k + 1. The following 
statements hold 

(i) If0<t<k then d (Atom 1 (s(k))) is of the form ((0,2)1, l fc -*) if t is even 
and of the form (2, (0, 2)^2" , l fc ~*) if t is odd. 

(ii) Ifk+1 < t < 2k+l then d (Atom 1 (s(k))) is of the form (0, (2,0) k+1 ~ i ^- , l t ~ k ) 
if t is odd and of the form ((2, 0) fe+1 ~ ' ; l*~ fe ) if t is even. 

Consequently, the height of Atom t (s(k)) is equal to k + 1 if t is odd and k if t is 
even. 
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Proof. We prove this lemma by induction on t. For t — we have Atom°(s(k)) = 
(k, k — 1, . . . , 2, 1) so that d (Atom°(s(k))) = (l k ) which corresponds to the state- 
ment. 

For t = 1 then by definition of Atom Procedure one sand grain is added to the 
first column. The new configuration is Atom 1 (s(k)) = (k + 1, k — 1, . . . , 2, 1) and 
d (Atom 1 (s(k))) = (2, l fe_1 ) which corresponds to the statement. We assume that 
the statement holds till step t. We prove that it also holds for the step t + We 
consider the following cases: 

Case 1: If t ^ k and t is odd. We have 

d(Atom t (s(k))) = (2,(0,2)^,1*-') if t < k 

and 

d(Atom t {s(k))) = ((2,0) fc+1 -M*- fe ) if t > k + l. 
Since t + 1 is even, we does not add any grain at this step and just apply the PSPM 
rule. Hence, 

d(Atom t+1 (s(k))) = ((0,2)^, l fe -*-!) i{ t <k 

and 

d(Atom t+1 (s(k))) = (0,(2,0) fe+1 - i * i ,l t - fe+1 ) if t > k + l, 
respectively. 

Case 2: If i ^ k and t is even. By the definition of Atom Procedure one grain 
is added on the first column at the step t+1. It transforms (0,2)2, l fe ~* (resp. 
(2,(0,2)^,1^*)) into (0,(2,0) fc+1 -^ i ,l t - fe ) (resp. ((2, 0) fc+1 -3 , l*- fc )). 

Case 3: If t = k. We have then 

d (AW(s(fc))) = ((0, 2) k ) if t is odd 

and 

d(Atom t (s(k))) = (0,(2,0)^) if Ms even. 

We also have 

(Atom k+1 (s(k))) = (2,(0,2)^,1) and 

(Atom k+1 (s(k))) = ((2,0)^,1), respectively. 

So the forms of the associated configurations are described as the statement. 

Now it is remarkable that if the associated sequence of a configuration has the 
form ((0,2)', 1 TO ) then the configuration is of height 21 + m. So the rest of the 
lemma is straightforward. □ 

Corollary 2. The Atom Procedure transforms the configuration (k, k — 1, 2, 1) 
into the configuration (k + 1, fc, 2, 1) after 2k + 1 transitions. 

Next, we introduce two other procedures which are implemented on unimodal 
sequences 

Definition 7 (Alternating Procedure). Let a be a unimodal sequence whose highest 
column is placed at position and this column has a large enough number of grains 
to enable to distribute to its neighbors. An Alternating Procedure is a sequence of 
transitions by PSSPM rule such that 

(i) Column collapses on the right at odd steps. 

(ii) Column collapses on the left at even steps. 

Remark: 
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(1) For a given configuration, the Alternating Procedure will be not imple- 
mented forever. It stops when the column has not enough grains to 
collapse either on the left at some even step or on the right at some odd 
step. 

(2) At all steps in Alternating Procedure, all columns except column have 
only one direction to collapse due to the unimodality. Furthermore, at each 
step, the column choose exactly one direction to collapse. Therefore, 
the configuration obtained from a after t steps of Alternating Procedure is 
uniquely determined. 

(3) The column decreases exactly 1 after one step of Alternating Procedure. 
Consequently, after t steps of this procedure the column decreases exactly 
t grains. 



R 



R 



L 



Figure 7. 6 first steps of Alternating Procedure from (9). The 
arrow together with the direction R or L (Right or Left) corre- 
sponding to the direction along which the column (dark column) 
collapses 



Figure ^ illustrates that we can not implement more than 6 steps of Alternating 
Procedure from (9). 

Denote by Alt* (a) the configuration obtained from a after t steps of Alternating 
Procedure. The following facts are straightforward 

(Alt t (a)) >0 = Atom t (o >0 ) , 

and 

((^ t ( a )) <0 )- 1 = J 4tom t - 1 ((a' <0 )- 1 ), 

where (a' <0 ) 1 ^ s obtained from (<z<o) 1 by applying one step by PSPM rule. 
In other word, applying t steps of the Alternating Procedure on a is the same as 
applying t — 1 steps of Atom Procedure on a' <0 and t steps of Atom Procedure on 
a >Q . 

Particularly, if we start from the singe-column configuration then after some 
steps of Alternating Procedure the weights of the strict left part and strict right 
part by differ at most 1. 

We have the following lemma 
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Lemma 2. Let n, k be positive integers such that n > k 2 and let a — (1, 2, k — 
l,n — k 2 + k,k — 1, 2, 1) be a unimodal sequence. Then a is reachable from (n) 
after k 2 — k steps of Alternating Procedure. 

Proof. We prove this statement by induction on k. For k = 1 then the configuration 
(n) is reachable from itself. We assume that the configuration (1, 2, . . . , t — 1, n — 
t 2 +t, t — 1, . . . , 2, 1) is reachable from (n) after t 2 — t steps of Alternating Procedure. 
That means Alt^-^n) = (1, 2, . . . , t- 1, n-t 2 + t, t- 1, . . . , 2, 1). We need to show 
that the unimodal sequence (1,2, ... ,t,n — t 2 — t,t, ... ,2,1) is reachable from (n) 
by Alternating Procedure (where n > (t + l) 2 ). 
Since 

(Alt t2 -\(n))) <0 = (^/-*((n))) >0 = s(t 1), 

by Corollary || and the above remark it is sufficient to show that the Alternating 
Procedure can be implemented on Alt* (a) in 2t steps. 
Thus, for 1 < i < 2t then by Lemma [j] we have 

h(Atom i ^ 1 (Alt t2 -\a)) <Q ) <t, 

and 

h{Ato7n l ~ 1 {Alt t2 ~ t (a)) >a ) <t 

In addition, due to n > (t + l) 2 we have 

n - {t 2 - t) - (i - 1) > n - (t 2 - t) - {2t - 1) = n - i 2 - t + 1 > t + 2. 

So after i — 1 steps of Alternating Procedure the column of Alt t2 ~ t (a) is still 
collapsible. This completes the proof. 

□ 

Corollary 3. Let k be positive integer. Then the unimodal sequence (1, 2, . . . , k — 
1, (k), k — 1, ... ,2,1) is reachable from (k 2 ) in PSSPM by Alternating Procedure. 

Let s(n, k) = (1, 2, . . . , k — 2, k — 1, (n — k 2 ), k, k — 1, ... ,2,1) be a unimodal 
sequence whose the column is of height n — k 2 , the left strict part by is a stair of 
k — 1 steps and the right strict part by is a stair of k steps. We have the following 
lemma 

Lemma 3. Let n, k be positive integers such that n > (k + l) 2 + (k + 1). Then the 
Alternating Procedure transforms the configuration s(n, k) into the configuration 
s(n, k + 1) after 2k + 1 steps. 

Proof. By Lemma [l], it is sufficient to prove that the column of s(n, k) has a large 
enough number of grains to enable to collapse in 2k + 1 steps. We prove this by 
induction on the number of implemented steps. For 1 < t < 2k we have 

h(Atom t (s(n,k) >0 )) < k + 1 and h(Atom t ((s(n,k) <0 y 1 )) < k + 1. 

On the other hand, since n > (k + l) 2 + (k + 1) we have 

n - k 2 - (t - 1) > n - k 2 - (2k - 1) 
>(k+l) + 2 

> h(Atom t (s(n, fc)>o)) + 2, 
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and 

n - k 2 - (t - 1) > (k + 1) + 2 

> h i^Atom 1 ((s(n, fc)<o) _1 )) + 2. 

It implies that the column of s(n, k) can collapse at the step t. That means the 
Alternating Procedure can implement in 2k steps. 
Moreover, by Lemma [l] we have 

(Alt 2k {s{n, k))) >Q = Atom 2k {s{n, fc)>o) = (k,k, . . . ,2, 1) 

and 

((Alt 2k (s(n, fc))) <0 ) ~* = Atom 2 "- 1 ((s(n, fc)< )" 1 ) = (fe, fe — 1, — ,2, 1), 

(since (s(n, fc) <0 ) 1 is stable, after one step by PSPM rule on (s(n, fc)<o) _1 we 
obtain itself). 
Besides, 

h (Alt 2k (s{n, k))) = (n-k 2 )-2k>k + 2 = h ((Alt 2k {s(n, fc))) >0 ) + 2. 

So that we can apply the Alternating Procedure on s(n, k) in 2k + 1 steps and 
obtain s(n, k + 1). □ 

Now we define the second procedure implemented on unimodal sequences. 

Definition 8 (Pseudo- Alternating Procedure). Let t be a positive integer and let a 
be a unimodal sequence. We assume that the highest column of a is placed at position 
and this column has a large enough number of grains to enable to distribute to 
its neighbors. The Pseudo- Alternating Procedure of t steps on a is a sequence of 
t transitions by PSSPM rule starting from a such that Alternating Procedure is 
applied from step i 2 + 1 to step (i + l) 2 (for i — 0, 1, ... , — 1) and from step 
[l/i] 2 + 1 to step t. 




Figure 8. 9 first steps of Pseudo- Alternating Procedure on (13) 
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Remark that a Pseudo- Alternating Procedure is a concatenation of Alternating 
Procedure, however it is not an Alternating Procedure (see Figure ||). 

Denote by PAlt 1 (a) the configuration obtained from a after t steps of Pseudo- 
Alternating Procedure. 

Figure |] shows 9 steps of Pseudo- Alternating Procedure starting from (13). The 
dashed area illustrates that Alternating Procedure is applied from the step 2 to 
step 4 corresponding to the case i = 1 in the above definition. It is also noticeable 
that we are not able to continue implementing Pseudo- Alternating Procedure on 
(13) more than 9 steps although after 9 first steps the darked column of height 4 is 
still collapsible on the left but it is un-collapsible on the right. 

Lemma 4. Let n,k be positive integers such that n > (k + l) 2 + (k + 1). The 
Pseudo- Alternating Procedure transforms the configuration (n) into the configura- 
tion (1,2, ...,k,n- (k + l) 2 , fc + 1, 2, 1) after (k + l) 2 steps. 

Proof. We prove this by induction on k. For k = 1 then the configuration (n— 1, 1) is 
reachable from in) , where the column collapses on the right by Pseudo- Alternating 
Procedure. We assume that the configuration (1, 2, . . . , t— 1, n — t , t, t — 1, . . . , 2, 1) 
is reachable from (n) after t 2 steps of Pseudo- Alternating Procedure (t < k). We 
can write PAlt^ = (1, 2, . . . , t — 1, n — t 2 , t, t — 1, . . . , 2, 1). We need to prove that 
the configuration a = (1, 2, . . . , t, n — (t + l) 2 , t + I, ... ,2,1) is also reachable from 
(n) if n > (t + l) 2 + (t + 1). Due to the determination of the Pseudo- Alternating 
Procedure it is equivalent to prove that a is obtained from PAlt* 2 ((n)) after 2t + 1 
steps of Alternating Procedure. This is actually obtained from Lemma |[ □ 

Now given a fixed-point form P = {j>\,pi, . . . ,pk) of SSPM, we put 

Div i (P) = \w(P <i )-w(P >i )\, 

and 

Div(P) = min {Div\P)}. 

l<i<k 

We call Div(P) the symmetric difference of P. A column i of P at which Div l (P) 
gets minimum is called the symmetric separator of P. 



FIGURE 9. Two fixed-point forms P = (1,2,3,4,5,5,4,3,3,2,1) 
and P' = (1, 2, 2, 3, 4, 5, 5, 4, 3, 2, 2, 1) with their separators marked 
by dark columns, and Div(P) = 2; Div(P') = 5 

Remark: A fixed-point form may have more than one symmetric separator. For 
instance, the form (1, 2, 2, 3, 4, 5, 5, 4, 3, 2, 2, 1) contains two symmetric separators: 
column 6 and 7 (see Figure ^). 
We have the following lemma 

Lemma 5. Let P = (j)\,pi, ■ ■ ■ ,Pk) be a fixed-point form of SSPM and h — h(P). 
Then there exists an index i of P such that pi = h and Div l (P) < h. Moreover, 
for this position i we have Div % (P) = Div(P). 
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Proof. Put t = min{i : pi = ft}. Since P is a fixed point of SSPM, we have 



and 



ft(ft-l) , , h(h-l) , 



ft(ft-l) , x ft(ft-l) 

j <w(P >t ) < \ ' + 3ft. 



2 - v ~ 2 
So that < \P<t — P>t\ < 3ft. We consider the following cases 

Case 1. If < \P<t — P>t\ < h then t is the position satisfying the statement. 
Case 2. Ifh + 1< \P>t — P<t\ < 3/i, we claim that < + 1 is the position satisfying 
the statement. First we show that pt+i = h. On the contrary pt+i < h — 1, we 
have 

h(h-l) , . h(h-l) , 
1 2 ; < ^(^>t) < - 2 ^ + h - 1. 

Therefore, |P>j — P<t | < ft — 1 which contradicts the condition |P <t — P >t \ > ft + 1. 
So that we must have pt+i = h and this also implies that 

h(h - 1) , 
P>t> V 2 ; + fe>P< t - 

On the other hand, we have -P<(t+i) = -P<t + ft and P>(t+i) = -P>t — ft- Hence 

\P>(t+i) ~ P«t+i)\ = \P>t - P<t - 2h\ < ft. 

Last, we prove that Div l (P) = Div(P) by showing that for other positions t then 
Div t (P) > ft. Thus, it is sufficient to show this is true for positions t satisfying 
that pt < ft — 1. Since the strict part of P (which contains the column of height ft) 
by t is of weight at least h ( h + 1 1 corresponding to the sequence (ft, ft — 1, . . . , I) and 

the rest strict part by t is of weight at most corresponding to the sequence 

(ft - 1, . . . , 2, 1). We have 

t , , h(h+l) ft(ft-l) 

DivHP) > y ' i > ft. 

y ' ~ 2 2 ~ 

This completes the proof. 

□ 

Now we present the proof of Theorem || 

Proof in details of Theorem |^. We recall here some notations used in the main idea 
of the proof Theorem ||: 

d = Div(P) and ft = ft(P). 

Now we represent precisely the construction the sequence of transitions by PSSPM 
rule to obtain P from (n) placed at 0. 

(1) Applying the Pseudo- Alternating Procedure from the step to the step d 2 . 

(2) Applying the Alternating Procedure from the step d 2 + 1 to the step n — ft. 

(3) Applying the deterministic procedure from the step n — ft + I (from this 
step there is no any column collapsible on both sides). 

It is noticeable that the construction showed above also works for both P symmetric 
and not symmetric. We recall that we are not always able to apply the Pseudo- 
Alternating Procedure and the Alternating Procedure. The next is to prove that 
we can implement d 2 steps of Pseudo- Alternating Procedure on (n); then we can 
implement (n — ft — d 2 ) steps of Alternating Procedure on PAlt d (n); last, we get P 
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at the end of the evolution. These statements will be proved in the following three 
lemmas. 

Lemma 6. It is possible to implement d 2 steps of Pseudo- Alternating on (n). 

Proof. By Lemma ||, it is sufficient to show that n > d 2 + d. Since P is a fixed point 
of height h and of symmetric difference d, so P must contain at least h 2 + d grains 
corresponding the configuration (1, 2, . . . , h— 1, h, h— 1, . . . , d+1, d,d, . . . ,2, 1). By 
Lemma |^, we have d < h and so n > d 2 + d and 

PAlt d ' 2 ((n)) = s(n,d). 

□ 

Lemma 7. It is possible to implement (n — h — d 2 ) steps of Alternating Procedure 
on PAlt d2 ((n)). 

Proof. On the contrary we assume that t is the first step at which this does not hold 
and 1 < t < n — d 2 — h. By Lemma || and the definition of Alternating Procedure, 
we have the following facts 

(Alt* (s(n,k))) Q = n-d 2 -i, 

(A/t < (s(n,fe)))_ 1 < h + 1 and (Alt i (s(n,k))) 1 <h + 1, 

for all i = 1,2, . . . ,t - 1. 

Put Alt t ~ 1 (s(n, d)) — b. Since we can not implement the Alternating Procedure 
on PAlt d ((n)) at the step t, we must have 

fro < b-i + 1 < h + 2 in case t is even 

and 

bo < bi + 1 < h + 2 in case t is odd. 

Combining with the hypothesis of the contrary t < n . — d 2 — h, we deduce for both 
cases that 

h+l<n~d 2 ~t+l = b <h + 2. 

Hence, either n — d 2 — t+l = h + lorn — d 2 — t + l = h + 2. 

On the other hand, we have w(P > o) — w(P < o) = d and w(P < o) + w(P > o) + h = n, 
so that n — d — h is even. Thereby, n — d 2 — t + 1 = h + 1 if t is even and 
n — d 2 — t+\ = h + 2\it\s odd. We consider these two cases 

Case 1. t is even. We have b = h+1 and 6_i > h. Since there is no grain added 
on the left of Alt l ~ 2 (s(n, d)) at the step t — 1, we get 

(Alt t - 2 (s(n,d)))_ 1 > h. 

It is remarkable that in the Alternating Procedure, the strict left and the strict right 
by sequentially fulfill the stairs before creating the new stairs of greater length. 
So that as from the step d 2 + 1 — 2 the strict left part is at least the stair (1, 2, . . . , h) . 
This is also true for the strict right part by since its weight is more than that of 
the strict left part. So during the evolutions of applying the transitions arbitrarily 
by the PSSPM rule on 6, the column of 6 never collapses again. Therefore, at 
the end of these evolutions, we get fixed points having the same height h + 1. We 
will show this is a contradiction. We assume that P' is a such fixed point. Then 
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(ft + l)(ft + 2) 



w (P>o) ~ w (P<o) = d+1, and by Theorem [| P and P' are of representation as 
follows 

P' = (l, 2, . . . , a', a' . . . , ft + 1, . . . , a' + d + 1, a' + d + 1, . . . , 2, l) 

and 

P = (l,2, ... ,a,a ... ,h, ... ,a + d,a + d, ... ,2,1). 

Since w{P') = w(P) and P ' = P + 1, it implies that a' = a - 1 and Pi j = P( = ft. 
So that Pi = P-x = Pq = h and P contains at least 3 different plateaus of 4 
following plateaus: P-h+a-2P-h+a-i of height a, P-iPq of height ft, PqP\ of height 
ft and Ph- a -d+\Ph-a-d+2 of height a + d (and ft — a — d > 0). This contradicts 
the condition that P is the form of a fixed point of SSPM . 
Case 2. t is odd. We have 

b = h + 2 and b\ = h + 1 . 

Since there is no grain added on the right of j4ii t_2 (s(n, d)) at step t — 1, we have 

(Alt t - 2 {s(n,d))) 1 >h+l. 

Therefore, 

w)(^ t - 2 (s(n,d)) J>0 ) > 
and 

w(P) = w(A/t^ 2 (s(n, d))) >(h+ l)(h + 2)-d + h + 2 

>h 2 + 4:h + 4-d>h 2 + Sh. 

This contradicts the condition that P is one fixed point of SSPM(n) and of height 
h (since h 2 + ih corresponds to the configuration (1, 2, . . . , h — 1, ft, h, h,h,h — 
1,...,2,1)). ' □ 

In summary, we conclude that starting from (n) we can implement the Pseudo- 
Alternating Procedure in d 2 steps then the Alternating Procedure in n — d 2 — h 
steps. As a result we obtain the unimodal sequence whose all columns are of height 
less than or equal to ft. 

Lemma 8. The procedure applying the transitions by PSSPM rule on the config- 
uration obtained after the processes (1) and (2) above is deterministic. Moreover, 
it converges to P. 

Proof. This is straightforward from the fact that all columns of the configuration 
obtained after two above processes are of the height less than ft. Furthermore, 
the column of the configurations in the procedure is always of the height ft. 
Therefore, there is no any columns collapsible on both sides and the procedure is 
deterministic. On the other hand, the left side part and the right side part by 
evolve independently by PSPM (so it is deterministic) to reach their unique fixed 
points. The difference of two these parts always is d. So at the end of the procedure 
we get P. □ 

The three above lemmas have ended our proof of the theorem. □ 

Figure [To| shows one way to obtain the fixed point (1, 1, 2, 3, (4), 3, 3, 2, 1) (of form 
112343321) from the single column configuration of 20 grains. The dashed zone 
illustrates 4 steps of Pseudo-Alternating Procedure to obtain the s(20, 2) (which 
is 1(16)21). The dotted zone illustrates 12 next steps of Alternating Procedure to 
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obtain the configuration of height 4 at column (equal to the height of the given 
fixed point). Last, it needs to apply 3 steps of PSSPM rule when there is no 
column collapsible on both sides to obtain the final fixed point. 

We also remarkable that the evaluations of inequalities in Lemmas ^, |!||| and 
specially Lemma ^ are very sensible. Furthermore, our procedures we constructed 
above are not commutative each other in general, especially when the center column 
has not enough grains to distribute to its neighbors. So this may be easy to lead 
to another fixed-point form of SSPM not the one we expect. For instance, for 
the form P = (122221) we have h(P) = 2, d(P) = 2 and if we start from (10) 
and implement the way shown in proof of Theorem ^| then we get exactly P. But 
if we first do the Alternating Procedure in 4 steps (which is equal to n — d 2 — h) 
and next do the Pseudo- Alternating Procedure, then we are able to implement this 
procedure only in 3 steps not 4 steps (which is equal to d 2 ). Hence, we final get 
the form 123211 which is not 122221. 



\ (20) R -+. (19)1 (18)2 _ L _^ 1(17)11 _ R 1(16)21 



R , x L 
_ . 1(15)31 . 



2(14)22 3+ . 123(5)4221 _ 124(4)3321 133(4)3321 



_ ^ 223(4)3321 _ _ ^ 1123(4)3321 

Figure 10. A way of length 
(1,1, 2, 3, (4), 3, 3, 2,1) from (20) 
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Last we give a upper bounded for the shortest length of the ways to reach a fixed 
point in PSSPM. 

Corollary 4. Let Tpsspm(tt-) be the shortest length to reach a fixed point in 
PSSPM(n). Then n — [y/n] < TpssPM{n) < n. 

Proof. Recall that if h is the height of a fixed point of PSSPM{n) then h = [y/n\ 
or [\/n\ — 1. Hence, TpsspMlji) > n — h. From the way we constructed above, it 
takes n — h transitions of applying Pseudo- Alternating Procedure and Alternating 
Procedure; then it takes at most h transitions of applying the final procedure to 
reach a fixed point of PSSPM(n). Therefore, TpssPMM — n - ^ 



4. Conclusion and Perspective 

We proved that beginning with a singleton column of sand grains, the sequential 
model and the parallel model produce the same fixed-point forms. To tackle the 
problem, for each fixed-point form of SSPM, we construct an explicit way of 
PSSPM evolution to obtain this fixed-point form. Every configuration in this 
way has a "smooth" form, even it can be characterized by a formula on the time 
of the evolution, whereas it is difficult to capture the forms of general reachable 
configurations. 
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Actually, the problem of finding a shortest way to reach a given fixed-point form 
of SSPM is interesting to explore. The way we constructed is not always a shortest 
way although it reveals many interesting properties to be possibly a shortest way. 
In fact, the difference between the length of our constructed way and the one of the 
shortest ways is at most [y/n\- We do not know so far an explicit formula of the 
length or the behavior of such shortest ways. 
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Figure 11. The right-furthest fixed point with n = 144 

During the PSSPM evolution we constructed, the original column is always a 
highest column, and it never receives any grain from its neighborhoods. It would be 
interesting to investigate the problem in which the positions of the fixed points are 
considered. In this problem, the fixed points of PSSPM are not the same as those 
of SSPM. All the fixed points of PSSPM might be fully characterized by the 
furthest fixed points (the maximum and minimum fixed points with respect to the 
lexicographic order). A possible way to obtain the right- furthest fixed point is that 
at a current configuration, each column always collapses on the right if it is possible. 
By doing the experiments on computer, it is surprising that when n — (8k + 4) 2 
for some k 6 N, the furthest fixed point has a nice pyramid-shape which has no 
plateau and the right-most grain is at distance Ilk + 4 . For example, let n = 144. 
Then the furthest fixed point is illustrated Figure |ll|. It is reasonable to come up 
with the following conjecture 

Conjecture. ~ where d(n) denotes the distance of the right-most grain 

in the right-furthest fixed points to the original column. 

Up to now, we do not know how to prove the above conjecture or disprove it. 
Maybe, it needs a deeper analysis on the whole space of reachable configurations 
of PSSPM. 
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